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ABSTRACT 

Two different model-independent mapping techniques have been applied to CGRO 
OSSE, SMM, TGRS and balloon data and reveal a feature in the 0.511 MeV e+-e" 
annihilation radiation pattern of our galaxy centered at £ ~ —2° and h ~ 10° with 
a flux of ~ 5 • 10~^ 0.511 MeV ph cm~'^ s~^. If near the galactic center, then e"*" 
sources are producing ~ lO^'^ e"*" s^^ which annihilate ~ 1-2 kpc above the galactic 
plane. A starburst episode within the inner few hundred pc of our galaxy would 
drive hot pair-laden gas into the halo, with the one-sidedness pointing to the site of 
initial pressure release at the onset of the starburst activity. Positrons lose energy 
and annihilate as they are convected upward with the gas flow, and we calculate 
high-latitude annihilation patterns and fluxes in accord with the observations. Changes 
in the ionization state when the escaping gas cools could give annihilation radiation 
substructure. The fountain of hot (~ 10^-10^ K) gas rising into the galactic halo would 
be seen through its enhanced dispersion measure, thermal emission, and recombination 
radiation. 

Subject headings: Gamma Rays: Theory — ISM: Supernova Remnants, Jets and 
Outflows — Nuclear Reactions, Nucleosynthesis, Abundances 
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1. Introduction 

In addition to the previously identified (Purcell et al. 1993; Ramaty, Skibo & Lingenfelter 

1994) 0.511 MeV annihilation glow from the disk and nuclear bulge of our galaxy, recent analyses 
(Purcell et al. 1997a, 1997b; Cheng et al. 1997) of data obtained with the Oriented Scintillation 

Spectrometer Experiment (OSSE) on the Compton Gamma Ray Observatory, the Solar Maximum 
Mission, the Transient Gamma Ray Spectrometer (Teegarden et al. 1996) on the Wind spacecraft, 
and GRIS, FIGARO, and HEXAGONE balloon experiments show a very significant annihilation 
emission component ~ 5 10° north of the galactic plane in the general direction defined by the 
axis of the galactic center lobe (Pohl, Reich & Schlickeiser 1992). Imaging limitations of the 7-ray 
telescopes make it impossible to resolve fine structure, but the high latitude component appears 
extended rather than pointlike. If near the galactic center at distance Sdg kpc, then e"*" sources are 
producing ~ (2 — 3//2)^^4 • 10"^^ (i| e+ which annihilate ~ 1-2 kpc above the galactic plane. 
Here / is the fraction of e"*" which annihilate via positronium (Ps) formation. 

Positron sources include radioactive emitters from supernovae and novae, black hole jets, 
and low-energy (~ 10-100 MeV nuc^^) cosmic rays. Due to proximity to the galactic center and 
the feebleness of annihilation radiation from star- forming regions, we conclude that the region of 
enhanced e+-e~ annihilation radiation exists about 1 kpc above the galactic center. A starburst 
episode within the inner few hundred pc of our galaxy would drive hot pair-laden gas into the halo, 
with the onc-sidedness pointing to the site of initial pressure release at the onset of the starburst 
activity. Here we investigate the production, transport, and annihilation of e"*" convected to high 
galactic latitudes in a wind produced in the galactic center region. 

In Section 2, we consider e"*" production from supernovae, and summarize evidence for 
starburst activity in the galactic center region. A model for e+ transport and annihilation in a 
fountain model is described in Section 3, and used to calculate a map of the annihilation emissivity 
from the direction of the galactic center. Predictions and a summary of the results are given in 
Section 4. 

2. Positron Production from Supernovae and Galactic Center Activity 

The origin of positrons through the decay of radioactive nuclei is confirmed by observations 
of ^^Co and ^^Co nuclear decay lines from SN 1987A (Matz et al. 1988; Kurfess et al. 1992), 
the ^^Al 1.809 MeV line found in clumped structure along the galactic plane (Diehl et al. 

1995) , and the ^^Ti nuclear decay line from Gas A (lyudin et al. 1994). The mean lifetimes of 
the ^^M^ssqq j^^j 56co^56pg decays in the ^^Ni-^^eCo-^^^Fe chain are 8.8 and 111.4 days, 
respectively, with a /3+ emitted 19% of the time in the latter reaction. Type la supernovae 
involving white dwarf detonantion or deflagration produce a time-averaged e"^ production rate 
jyS^sja ^ ^ 10^^?7_2M56-/Via/c s^^, whcrc M56 is the average number of Solar masses of 
synthesized ^^Fe per SN la, Nia,/c is the number of SN la per century throughout the Milky 
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Way, and rj = W~'^7]^2 is the escape fraction. Chan & Lingenfelter (1993) calculate the 
escape fraction to be in the range 0.1 ^ r/_2 ^ 10, and report (Woosley & Weaver 1992) that 
0.6 < M56 < 0.9. For core-collapse Type II supernovae, only an average O.O8M0 of ^^Fe is 
synthesized per SN, and the escape fraction could reach 0.7% for well-mixed ejecta, implying that 

^56,11 ^ 1042^_2(M56 /0.08)iVjj/C 6+ S'^. 

The mean lifetime of the ^'^Ti— >^^Sc— >'^^Ca chain is 78 yrs, with a positron produced 95% of 
the time in the latter decay. Compared to the mass of ^^Fe, a ^^Ti mass fraction = 10~^C_4 
in the range 0.3 ^ C-4 ^ 1-4 is calculated for SN la and C_4 < 25 is found for SN II (Chan &; 
Lingenfelter 1993). This gives N^l^ = 8 x 10^-^C-4-^56-^SN/Ci noting that most e"*" escape from the 
SN ejecta and mix with the surrounding medium. ^^Al decays into ^^Mg with a mean lifetime of 
10^ yrs, producing a positron 82% of the time. Production of ^^Al is most important in Type 
II SNe, though only 0.3 • 10"^ ^ Mae ^ 20 • 10"^ Solar masses of "^^Al are produced per SN II, 
depending on initial stellar mass and the treatment of semi-convection (Prantzos 1996). This gives 
jy2^6,ii ^ l0*^^(M26/10^^)iVji/c s^"*^, a value probably insufficient to account for the 1.5 M© of 
26 Al required to explain the total observed (11) 1.809 MeV ^'^Al line flux of 3.1(±0.9) • 10"^ cm-^ 
s~^ for which, at least, Wolf-Rayet stars and novae make important contributions. 

An episode of starburst activity in a region a few hundred pc across enclosing the galactic 
center could account for Ginga observations (Koyama et al. 1989; Yamauchi et al. 1990) of 6.7 keV 
emission from He-like Fe at temperatures ^ 10^ K. The hot gas in this scenario is produced by 10^ 
SNe over the preceeding 10^ yrs, which would be mostly SNe II from massive star evolution unless 
the activity had persisted for a period much exceeding 10^ yrs. Other indications (Hartmann 1995) 
of explosive events near the galactic center include large-scale X-ray structures, radio structures, 
and indeed the 1.809 MeV line. The radio emission of the galactic center lobe is explained by the 
synchrotron emission of nonthermal electrons convecting and diffusing outward (Pohl et al. 1992). 

Within the limits of uncertainty, (3~^ production from SNe II could account for the flux 
of the high-latitude annihilation glow if Ng^/Q ~ 1 in the central 100-200 pc nuclear region 
of our galaxy and, furthermore, if e"*" are transported to the annihilation site and efficiently 
annihilated. Pair-laden hot gas from the galactic center starburst would vent into the galactic 
halo while expanding, radiating and cooling, and slowing in transit from the galactic center to 
the lower-pressure galactic halo. While convecting outward, e"*" would cool and annihilate with 
the plasma electrons in the hot wind, and this is our explanation for the observations (Purcell et 
al. 1997a, 1997b; Cheng et al. 1997). The venting would be preferentially one-sided due to the 
location of the initial starburst activity and the morphology of the confining gas.[| 



^We note that either Sgr A*, the putative 10^ M0 black hole at the Galactic Center, or the less massive Einstein 
source (IE 1740.7-2428), weighing in at ~ 10-10'^ M0, could produce a large enough production rate if annihilation 
were efficient, but the origin of the elevated target region requires additional explanation. 
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3. Annihilation Fountain Model 

We sketch a model for the fountain of rising annihilating gas energized by a starburst 
episode.0 Suppose that gas rises with speed 10'' cm s~^ from a region of radius = lOOrioo 
pc undergoing a starburst phase, implying a gas crossing time of order rioo/i'? million years. We 
approximate the shape of the volume of the gas expanding into the galaxy's halo by an inverted 
cone with opening angle X) cut off at both ends, so that the cross sectional area of the fountain 
is •7rr^(l + u)^, where u = ztanx/rt and z is the height above the galactic plane. Our major 
simplification, which can be relaxed in more general treatments, is that the gas rises with constant 
velocity. Continuity of the mass flux for a steady-state situation considered here (see Ramaty 
et al. 1992 for a treatment of time-dependent injection and annihilation in a uniform medium) 
implies that the density of the gas at the base of the fountain is = O.13M0/c/(f7r^QQ) cm~^, 
where Mqjq Solar masses of gas are expelled per century from the starburst region and rise into 
the halo. The time- and spatially- averaged density distribution as a function of height z above 
the galactic plane is therefore np{z) = '^^/(l -|- u)^. 

The time scale for a positron injected with kinetic energy meC^(7 — 1) to thermalize with the 
background hot thermal gas is controlled primarily by Coulomb losses at mildly relativistic and 
nonrelativistic energies. Coulomb losses operate on a time scale of ~ 0.5/3(7 ~ l)/(^-i^3o) million 
years, where A30 is the Coulomb logarithm divided by 30, the density is 0.1n_i protons cm~^, and 
/3c is the positron's speed. The kinetic energy distributions of positrons entering the ISM after 
SN II explosions are given for the ^^Co— >^^Fe and the ^'^Sc— >^^Ca decays by functions peaking 
near 0.6 MeV with FWHM widths of ^ 0.9 MeV and high-energy tails reaching to « 1.45 MeV 
(Chan & Lingenfelter 1993). The fraction of e"*" which annihilate in flight prior to thermalizing 
usually amounts to < 10% and these positrons do not contribute to the 0.511 line emission 
(Murphy, Dermer & Ramaty 1987). After entering the thermal pool, the positrons annihilate in a 
fully ionized thermal plasma through direct annihilation and radiative combination on a thermal 
annihilation time scale of 30/(A_i4n_i) million years, where the annihilation rate coefficient 
A = 10~^^A_i4 cm^ s~^ (Bussard, Ramaty & Drachman 1979). When a significant fraction of 
neutral atoms or partially ionized ions are present, the annihilation rate can dramatically increase 
because A_i4 — > 10^ near temperatures of ~ 10^ K due to the onset of charge exchange processes 
which have atomic-sized cross sections. A pair of 0.511 MeV line 7 rays comes from direct 
annihilation in the thermal gas and, one-quarter of the time, from annihilation via Ps formation. 

We calculate the annihilation flux at height z after tracking the energy evolution and 
spatial propagation of the e"*" following injection. The e"*" height-dependent injection function is 
approximated by a Gaussian function which is allowed to be, in general, offset from the galactic 
plane. The FWHM width of the injection function is taken to be 180 pc, comparable to twice the 
scale height of massive stars. Though diffusion must be important for very relativistic e^ (Lerche 



^The general picture also applies to plasma ejection from black hole jets. 
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&; Schlickeiser 1980), we assume that the nonthermal and thermal e+ are entrained in the hot gas 

and convect away from the galactic midplane with constant speed vq. The general e"*" equation 
of motion involves Coulomb, bremsstrahlung, adiabatic expansion, synchrotron and Compton 
energy losses, but for injection, only Coulomb and adiabatic expansion losses are important. 
The e"*" energy-loss rate from adiabatic expansion is given by — 7adia = (47 — 3 — j~^)V /{3V), 
which bridges the nonrelativistic and relativistic regimes and applies to high-beta plasmas, which 
is suitable for magnetic fields weaker than ^ 10~^ Gauss. The volume expansion rate for the 
fountain geometry is V /V = 2vQtanx/[i"bi'^ + u)], from which we find that Coulomb losses 
dominate adiabatic expansion losses when /9(47 — 3 — 7"^) ^ K = 36A3o-MQ/c/b7'"ioo(tanx/0.1)]. 

After emission, positrons convect away from the galactic plane and, if injected with 
sufficiently low energies, thermalize with the background gas through Coulomb processes (excepting 
those few which annihilate in flight). We call the distance between injection and thermalization 
the Maxwell-Boltzmann length (MBL), which is a nonlinear function of the injection height Zi 
and the e+ injection kinetic energy meC^(7 — 1). When Coulomb processes dominate, the MBL is 
given by umb = [(1 + Ui)~'^ - K-^{(5i'yi - arccos7~^)]~^ - 1, where nMB(i) = -^mbO) tanx/rfo. The 
constant Kc = 2.4M0/cA3o/(f7r-iootanx). 

When K » 1, as with standard parameter values, positrons injected through (3'^ production 
thermalize close to their injection site. Thus the spatial dependence of positrons thermalizing 
with hot gas is, in this regime, equal to the energy-integrated /^"'"-injection function. Following 
thermalization, a positron continues to convect upward into the galactic halo until it either 
annihilates or merges with the dilute halo gas. The decay law through annihilation for the rising 
positrons is governed by the value of the temperature-dependent reaction rate coefficient A(r), 
the ionization state, composition, and density of the medium. For hot gases with T ^ 10^ K, the 
2-photon direct annihilation channel is most important with A_i4 = 1. The decay law for thermal 
positrons annihilating in a thermal gas is —N^{t) = N-^-{t)Xnp{t). From this equation, we derive 
the height-dependent differential production rate of 0.511 MeV annihilation line photons, given by 

2r)0 \(T) rz 1 1 

Here N^{z) is the MB injection function differential in height z, which is obtained by convolving 

the Z?"*" injection function with its MBL and integrating over the initial energies of the (3~^ positrons. 
The constant Ka = 0.039A_i4MQ/c/(?;7riootanx), nP, is the electron density at the base of the 
fountain, and u^'^ = z^'^ tanx/'^fo. 

In the approximation where e"*" thermalize close to their injection site, equation (1) is easily 
solved to give the results shown in Fig. 1. The total e"*" injection rate is lO^^A^+42 e"*" s~^ with 
iV+42 = 1. The solid curves give the integral 0.511 MeV photon production rate between the 
galactic midplane and height z, and the dashed curves give the differential 0.511 MeV production 
rate in units of 0.511 MeV ph s~^ z(pc)~^. The dotted curves represent the spatial 6+ distribution 
function, which are centered at and 50 pc above the galactic midplane in Figs, la and lb, 
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respectively. 

Before interpreting Fig. 1, note that iVo.511 MeV 

(< 00) ^ 2 ■ lO''^ 0.511 MeV ph s'^ if all e+ 
annihilate. This limiting value is reached only if most /3+ injection occurs high above the galactic 
plane so that few e"*" are convected to negative values of z. Standard parameter assignment with 
no source offset and tanx = 0.1 (Fig. la), corresponding to a 6° fountain opening angle, yields 
an integral annihilation flux of ~ 3 • 10^^ 0.511 MeV ph s~^, implying a one-sided annihilation 
efficiency of ~ 15%. In contrast, when the fountain opens to 45° (tanx = 1), the annihilation 
efficiency plummets because of the severe thinning of gas density with height. In Fig. lb, the wind 
speed parameter for the hot gas is increased by an order of magnitude from the standard value 
and the injection Gaussian is offset by 50 pc. This also reduces the annihilation efficiency because 
flux continuity demands a reduced density in a steady-state approximation. 

The overall shape of the differential 0.511 MeV production function is unusual in that it 
exhibits a broad plateau between ~ 50 pc and 1-2 kpc. Rather than decreasing monotonically 
with height above the galactic midplane, the maximum of the differential 0.511 MeV annihilation 
flux peaks downstream from the maximum of the Gaussian injection function. This can be 
understood by superposing contributions from discrete sources. The annihilation flux decreases 
downstream from a steadily emitting e"*" source as the radioactive debris convect outward with 
constant velocity. But the divergence of the wind vector at the galactic midplane means that the 
e^ density at a height zq is proportional to the number of sources between the galactic midplane 
and zq. This quantity increases with height above the galactic midplane for reasonable injection 
distributions until the source injection function tails off. The underlying assumption of time- 
and spatial-averaging appears to be satisfied for our system, because ~ 10^ SNe occur during 
the gas-crossing time scale of the starburst region. More detailed calculations must consider the 
spatial distribution and temporal evolution of the discrete SN events exploding throughout the 
region. 

The rate at which 0.511 MeV annihilation photons are emitted between ~ 0.1 and 2 kpc 
is ^ 2 • lO^iiV+42 and ^ 3 • 0.511 McV ph s'^ for the top curves in Figs, la and lb, 

respectively. To match the OSSE/S'MM/TGRS/balloon data analyses (PurccU ct al. 1997a, 1997b; 
Cheng et al. 1997) implies that annihilation occurs at the rate of pa 4 ■ lO^^dg 0.511 MeV ph s~^. 
Thus injection rates of iV_|_42 ~ 20 and7V-|-42 ~ 13 are required to match the observations for these 
two cases. Such rates could be supplied by /J"*" production from freshly synthesized ^^Ti in SN II 
given the uncertainty of the ^"^Ti mass fraction, the SNe rate, and the parameters of the outflowing 
wind. This rate could perhaps also derive from (3^ production in the "''^Ni— >^^Co^^^Fe chain if 
either the calculated escape fraction or amount of synthesized ^^Ni per SN II were underestimated. 
Black hole e"'' production could provide an additional source of the positrons. 

Fig. 2 shows a contour map of the annihilation flux produced by the galactic fountain, added 
to a disk and galactic bulge distributions fitted to the measured (Purcell et al. 1997a) annihilation 
emissivity. The axis of the fountain is directed away from the center of our galaxy along the 
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direction of the galactic center lobe, which is inclined in projection by ~ 20° from the axis of the 
North Galactic Pole. (Note that the location of the centroid of the spheroidal bulge flux is offset 
from the location of the galactic center.) The agreement of this idealized model with the data is 
satisfactory, and predicts that peak enhancement of the fountain's annihilation flux occurs ^ 100 
pc above the galactic plane. 

4. Predictions and Summary 

Due in large part to gamma-ray observations with the OSSE instrument on Compton, a 
new component of the ISM has been discovered: hot plasma pushing into the galactic halo from 
a region of starbiirst activity surrounding the galactic center a few hundred pc in extent. The 
existence of the outflowing wind is suggested by radio and X-ray observations (Pohl et al. 1992; 
Koyama ct al. 1989; Yamauchi et al. 1990; Morris &: Serabyn 1996) and represents a low-power 
analog of nuclear starburst outflows detected in M82 and NGC 253 (Shopbell & Bland-Hawthorn 
1997). Hot ionized plasma can be detected through pulsar dispersion measurements, and we 
predict a dispersion measure jump of ~ 50 when scanning across the fountain several degrees above 
the galactic plane. The difficulty is to have independent distance measurements of pulsars located 
> 8 kpc away. Continuum free-free and recombination Lya lines of H, He, and Ps are predicted 
from the fountain, but obscuring foreground gas makes it difficult to view the galactic center 
region at optical and UV wavelengths. As ^^Al is convected upward with the flow, a characteristic 
decay law for the diffuse 1.809 MeV emission (Chen, Gehrels & Diehl 1995) should be observed, 
though at a flux level requiring at least INTEGRAL telescope sensitivities. Turbulent hot gas in 
the central region might also explain the broadened 1.809 MeV lines observed (Naya et al. 1996) 
with the GRIS instrument. The width of the 27 0.511 MeV line from annihilation in hot gas will 
be broader than the galactic disk 0.511 MeV line emission because of the larger temperature of 
the plasma, depending in detail on the dust content of the fountain, and the 87 Ps continuum 
fraction / will be spatially varying (Ramaty, private communication). 

If the annihilation fountain is our first clear view of a conduit of hot gas venting the starburst 
activity near the galactic center into the galactic halo, then it is likely that cosmic rays are 
efficiently transported into and diffuse throughout the halo along this route. The > 100 MeV 
EGRET (Hunter et al. 1997) and 408 MHz Effclsbcrg (Haslam et al. 1982) all-sky maps show 
enhanced emission north of the galactic plane which could be ascribed to this effect. Positrons 
which fail to annihilate in the fountain will diffuse throughout the halo and annihilate if they 
encounter the gaseous disk. This could explain the large scale height of the disk component of 
the diffuse annihilation radiation ^ 10-20° away from the galactic center (Ramaty et al. 1994). 
Long-lived radioactivity produced by SNe would also trace this channel. Small-scale variations 
of the 0.511 MeV map would indicate localized OB associations, changes in the ionization state 
of the gas, and sites of black hole e"*" injection. More observations and analyses are required to 
understand better the activity occurring near the center of the Milky Way. 
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Fig. 1. — Height distribution of the differential (dashed curves) and integral (solid curves) 0.511 
Me V annihilation emissivities produced in a fountain of hot gas rising upward with constant velocity 
Vf). Except where noted, the opening angle of the fountain is x = arctan(O.l) = 5.°7, vq = 10^ V7 
cm s~^ with v-j = 1, and the radius of the base of the fountain is 100 pc. We assume that 1 Solar 

m,ass of hot gas is ejected into the fountain per century. The dotted curves show the spatial 
injection function normalized to a total injection rate o/lO^^ s~^ . Fig. 1(a) compares the effects 
of different opening angles x for an injection function symmetric about the galactic midplane. Fig. 
1(b) compares the effect of varying gas speeds for an injection function offset from the galactic 
midplane by 50 pc. 
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Fig. 2. — Contour plot of the model annihilation emissivity of the Milky Way including the disk, 
galactic bulge, and fountain component. The contours are in units of 10-(2+")/5 0.511 MeV ph 
cm~^ s~^ sf~^, with the central contour corresponding to n = 1. 



